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1 Abstract 

The powerful and highly collimated jets observed in active galactic nuclei and fi— quasars 
are likely to be connected to the accretion phenomenon via disks. Based on theo- 
retical arguments and quasi-stationary radiative MHD calculations, a model for an 
accretion-powered jet is presented. It is argued that accretion disks around black 
holes consist of 1) a cold, Keplerian-rotating and weakly magnetized medium in the 
outer part, 2) a highly advective and turbulent-free plasma inside r tr = 10 — 20 
Schwarzschild radii, where magnetic fields are predominantly of large scale topology 
and in excess of thermal equipartition, and 3) an ion-dominated torus in the vicin- 
ity of the hole, where magnetic fields undergo a topological change into a monopole 
like-configuration. The action of magnetic fields interior to r tr is to initiate torsional 
Alfven waves that extract angular momentum from the disk-plasma and deposit it 
into the transition layer between the disk and the overlying corona, where the plasma 
is dissipative and tenuous. A significant fraction of the shear-generated toroidal mag- 
netic field reconnects in the transition layer, thereby heating the plasma up to the 
virial-temperature and forming a super-Keplerian rotating, and hence centrifugally 
accelerated outflow. 

The strong magnetic field in the transition layer forces the electrons to cool rapidly 
which, in combination with the fast outward-oriented motion, yields a two-temperature 
ion-dominated outflow. The toroidal magnetic field in the transition layer is in ther- 
mal equipartition with the ions, whereas the poloidal component is in equipartition 
with the electrons. Such a strong toroidal magnetic field is essential for increasing 
the jet-disk luminosity in the radio regime. These gravitationally unbound outflows 
serve as seeds, possibly, for all the powerful electron-proton jets observed in accreting 
systems containing black holes. 
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2 Introduction 



Based on observational data, most of the systems containing jets are considered to be 
accreting systems, most likely with an accretion disk surrounding the central object. 
The most powerful and highly collimated jets are found to emanate from accreting 
systems containing, probably, black holes, with a maximum power attainable if the 
central objects are Kerr black holes (-BHs) rotating at their maximum rate Be- 
side the numerous extragalactic radio sources with high gamma- factors (7 > 3), the 
underlying engines powering the jets in microquasars such as GRS 1915 + 105 and 
GRO 1655 - 40 are believed to be spinning BHs [2*7j . 

Although there has been a significant progress in recent years toward understand- 
ing the morphology, propagation and shock diagnostics of jet-plasmas, no consensus 
has been reached yet about their basic driving mechanisms. Nevertheless, there are 
at least three ingredients that appear to be necessary for initiating jets that gained 
theoretical and observational supports: 1) large scale magnetic fields (-MFs), 2) an 
accretion disk and 3) a central object which dominates the disk dynamics gravitation- 
ally. Magnetic fields in particular are considered to play a major role in the evolution 
of jets. The fact that many jets have comparable radio and bolometric luminosities 
hints to a significant fraction of the magnetic energy in their total power. 
In the first place MFs extract rotational energy from the underlying disk and deposit 
it into the plasma at higher latitudes, thereby initiating motions that are principally 
outward-oriented. 

At relatively large radii MFs are pre- dominant ly toroidal, and their main action is to 
force the rotating particles to collimate around the corresponding axis of symmetry. 
The central object must be relatively massive compared to the disk, so that its gravi- 
tational well is sufficiently deep for initiating motions with speeds comparable to the 
escape velocity near its surface. 

In this paper we present a model for initiating jets from accretion disks around BHs. 
The model is a modified version of the truncated disks advective tori model (TDAT; 
[To] , in which a special attention is given to the role of MFs. 

Our approach relies on theoretical arguments supported by the radiative two-temperature 
MHD numerical calculations. The model is based on the following assumptions: 1) 
The Balbus-Hawley instability acts as dynamo that amplifies the unordered MFs up 
to thermal equipartition, 2) The Parker- and BH-instabilities in combination with 
reconnection close a dynamo cycle, through which a large scale MF is generated, 
in the manner that Tout & Pringle (1992) suggested. Once the large scale magnetic 
field grows beyond thermal equipartition, the generation and dissipation of turbulence 
will be suppressed. 3) A significant fraction of the toroidal magnetic field (-TMF) 
lines re-connects in the transition layer (-TL): a process which is treated by adopting 
a turbulent diffusivity 77 mag , 4) Turbulent dissipation preferentially heats the heavy 
particles rather than the electrons, so that the two-temperature description for the 
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plasma can be applied [22]. Fig. 2 shows the main phases in the evolution of a 
accretion-induced jet. 

The paper runs as follows. In Sec. 3 we review several model for jet initiation, and 
outline the necessity for a new model. The structure of magnetized disk is discussed 
in Sec. 4. In sec. 5 the governing equations to be solved and the method of solution 
are described, while the corresponding results are presented in Sec. 6, and end up 
with the summary of the results in Sec. 7. 



Figure 1: A schematic distribution of the rotational velocity Q and density p across 
the disk with magnetic fields (MHD) and without (HD). The action of the equiparti- 
tion MFs is to transport angular momentum from the disk to higher latitudes, where 
the particles are centrifugal-accelerated outwards. 



3 The necessity for a new model 

Several models have been suggested for initiating jets[2T], 0, |3Tj - jSj, jSZI, |3Uj . 

|Hj; IE] M |2"2*| . [IH]- In the following we outline several properties of a few of 
these models: 

1. The Blandford& Payne model (1982, henceforth BP82) relies on magnetic ex- 
traction of angular momentum and rotational energy from an underlying cold 
and Keplerian disk, i.e., from a standard disk (-SSD, Extraction is most 

efficient if the poloidal magnetic field is of a large scale topology and appro- 
priately aligned to the disk-normal (i.e., the angle from disk-normal 8b must 
be larger than 30°). The second important result is that, at large radii, colli- 
mation is achieved through the action of the shear-generated toroidal MF. The 
difficulties of this model are: 
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(a) The disk is infinitely thin. This implies that Lorentz forces exerted on 
the disk-plasma would force the inflow to rotate sub-Keplerian, making it 
difficult for the plasma to be flung out by purely centrifugal means 38J. 
Furthermore, taking into account that toroidal magnetic field component 
Bt(z = 0) = and adopting the density profile p ~ r ~ 3 / 2 as BP82 sug- 
gested, any dynamo cycle will not succeed to amplify B p at z = to 
values beyond thermal equipartition. One possible way to obtain 1 f3 > 1 
is through a collapse of the central part of the disk while freezing its mag- 
netic flux. In this case 2 B p ~ r~ 2 , which implies that if B p does not change 
topology, accretion will be terminated and the whole inflow will turn into 
outflow with unacceptably large mass load. 

(b) It has been argued that even if the poloidal MFs intersects the disk with 
9b > 30°, thermal assistance is still required for the disk-plasma to over- 
come the potential difference parallel to the MF-lines. 

(c) There is no special treatment to the flow in the deep gravitational well of 
the central object, which is necessary for increasing the total energy per 
gram of the inflow, and convert it into outflows with high speeds. 

(d) The model relies on ideal MHD treatment to a cold Keplerian disk. How- 
ever, since (3 > 1, MFs will likely suppress the generation and dissipation 
of turbulence. Hence, the disk is likely to be thermal unstable. If the accre- 
tion rate is high, the disk becomes optically thick to Synchrotron/cyclotron 
radiation; gyrating electrons emit radiation on the cyclotron frequency, 
that in turn repeatedly upscattered by the hot electrons, establishing 
thereby a thermal equilibrium. Since this occurs on a time scale com- 
parable or even shorter than the dynamical one, and in the absence of 
other efficient heating sources, the disk undergoes a runaway cooling: the 
temperature decreases and reach the lower limit T min = hu pea ^/k, where 
^peak is the peak cyclotron frequency. This thermally induced collapse of 
the inner disk gains acceleration if the accretion rate is low. Here the 
low frequency photons emitted by the gyrating electrons escape from the 
system without being absorbed or scattered, and so no lower limit for the 
temperature can be constructed. 

2. The X-wind model (Shu et al. 1994). Here the magnetically driven-wind em- 
anate, neither from the disk nor from the central object, but from the region 
around the co- rotation radius, r cor , where the effective gravity vanishes. The 
model relies primarily on the poloidal MFs of the central object which is sur- 
rounded by an accretion disk, essentially without large scale magnetic fields. 

1 /3 = P ma g/-Pgas, where P mag = (B p + B^)/Att, see Eq. 8 for further details. 
2 $ = 2irr 2 B p « 27rr 2 S e = const. 
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This approach relaxes the winding up problem of the poloidal MFs, which has 
been encountered in the Ghosh & Lamb model (1979). Moreover, as outflows 
here originate from the region close to the surface of the star, their propagation 
speed is estimated to be of the order of the escape velocity. 
When extending this model to accretion flows around BHs, the following prob- 
lems arise: 

(a) It is not applicable to accretion flows onto Schwarzschild BHs, as these 
objects do not possess dynamically stable poloidal MFs. If such holes are 
viewed as extreme stars, then the X-wind model rule out indefinitely the 
possibility of jet initiation from weakly magnetized slowly rotating objects. 

(b) Balbus & Hawley instability amplifies weak MFs up to approximately ther- 
mal equipartition on the dynamical time scale. In the innermost region 
of the disk, MFs of the central object are beyond super-equipartition with 
respect to the thermal energy of the disk-plasma; generation and dissipa- 
tion of turbulence will be suppressed, the plasma in the X-region (see Fig. 
2b in PHI becomes cold, and the radial accretion will be terminated. 

(c) Almost all young stellar objects have been observed to rotate far below the 
beak-up velocity. Although the FU Orinois objects accrete at relatively 
high rates, they are considered to be slow rotators. Fast rotators, how- 
ever, show excessive magnetic activities, which is in line with the Parker 
instability [lOj. In the X-wind model however, MFs in the X-region are 
nearly in equipartition with the potential energy of the flow. Such strong 
magnetic tubes are likely to be gravitationally unbound, cannot be an- 
chored deep in the star and therefore they float upwards to the surface on 
the dynamical time scale. This occurs if they are generated in the convec- 
tion zone or in the overshoot layer between the convection zone and the 
underlying rigidly rotating core, similar to the solar dynamo 

3. Advection-dominated inflow outflow solutions -ADIOS (Blandford & Begelman 
1999), 

ADIOS are special case of ADAF solutions (advected dominated accretion flows) 
in which the accretion rate is allowed to adopt the self-similar profile M. ~ r p , 
where < p < 1. Thus, the accretion rate decreases inwards, thereby giving 
rise to a substantial outflow with large mass load. ADIOS model is another 
step closer than ADAF toward explaining the low luminosity AGNs. Indeed, 
there are several numerical calculations that confirm the inward decrease of the 
accretion rate (e.g., [IB])- However, these calculations did not rule out 

the possibility that the outflows obtained might be large scale circulation. On 
the other hand, in addition to the difficulties associated with ADAF solutions, 
ADIOS model: 
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Figure 2: A schematic illustration of the jet-disk connection. In the outer region 
of the disk MFs are weak and advected-around MFs by the fluid motions. BH- and 
Parker instability in combination with reconnection and inwards accretion amplified 
the MFs up to thermal equipartition. The MFs become of large scale topology and 
trun the disk at r tr into advection-dominated (part a). Interior to r tr (b), turbulent- 
generation and dissipation are suppressed, and torsional Alfven waves become the 
dominant angular momentum transporter (c). torsional Alfven waves carry angu- 
lar momentum from the disk and deposit it into the thin layer above, where ions 
rotate super-Keplerian, virial-heated by reconnection of the TMF-lines, and start 
to centrifugal-accelerate outwards (part d), The ions in the TL are associated with 
internal, rotational and TMF-energy to start collimation at r co \, where the magnetic- 
diffusivity becomes vanishingly small (part e). For clarity, a schematic description of 
the density profile versus radius of the model is shown. 
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(a) relies on one dimensional treatment to an intrinsically multi-dimensional 
phenomenon. The model does not provide answers to numerous questions 
related to jet-morphologies, such as the origin of the high speed of jet- 
propagation, the role of the central object, the magnetic field topology 
appropriate for initiating jets and etc. 

(b) It largely overestimates the mass load of jets. Basically, it predicts that 
almost all accreted matter through the disk at large radii will re-appear in 
the jet. When applied to the galactic center, an accretion rate that is about 
one order of magnitude larger than the upper limit would be necessary to 
obtain a partial fitting with the observed SED (Yuan et al 2002). Further- 
more, neither ADAF nor ADIOS models are able to reproduce reasonable 
SED-fitting with observations in the radio frequency regime. They seri- 
ously underestimate the radio luminosities in AGNs and black hole X-ray 
transients [TT]. [TU] . 

(c) The dynamical range of ADIOS is largely overestimated. To clarify this 
point, let r tr be the transition radius, at which an SSD makes a transition 
to ADAF (Fig. 2). If both solutions are thermally and dynamically stable, 
why should an SSD change into ADAF at large radii, or vise versa? 
Mathematically: in the absence of external forces and for given Dirich- 
lect boundary conditions, there is at most a single weak solution to the 
set of the corresponding linear systems of equations, which fulfils the en- 
tropy principle. Accordingly, at large radii the system of equations are 
nearly linear (i.e., they can be linearized as there are no mechanisms for 
generating shocks or strong gradients), and the conditions imposed at the 
outer boundary determine the solution uniquely. Note that the equations 
describing accretion flows in SSD or ADAF are predominantly Eulerian, 
except the angular momentum and energy equations that contain second 
order partial differential operators 3 . For example, if we set the angular ve- 
locity at the outer boundary to be sub-Keplerian and replacing Le by Ab 
(see Eq. 12), the SSD changes to ADAF. However, since there is no reason 
to expect an excessive source of heating at large radii, as in ADAF case, 
only cold standard disks are likely to survive there. These disks continue to 
maintain their thin geometrical structures all the way down to the central 
object, where other effects come into play, e.g., heat conduction, change of 
the MF-topology, noticeable change of rotation-law and manifestation of 
the magnetic field and the gravitational well of the central object. Indeed, 
Abramowicz et al. (1998) concluded that a transition from SSD to ADAF 

3 In spherical geometry they read Le = V • AfldVS, and Lg = V • ?7turVf2, respectively. E and 
Afld are the radiation density and the radiative diffusion coefficient as defined in Eq. 12, O and 
rjtur the angular frequency and coefficient of turbulent diffusion, respectively. 
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may occur in the vicinity of the BH. Hujeirat & Camenzind (2000a) found 
that the transition is thermally unstable, and the plasma in the outermost 
part of ADAF collapses dynamically, thereby extending the SSD down to 
several last-stable radii, where the disk truncates and an ion-dominated 
and highly edvective torus emerges. 

An additional evidence that rules out transitions at large radii are MFs. 
ADAF/ADIOS models assume that MFs are in equipartition with the ther- 
mal energy. If this is the case, then this should apply also when MFs are 
of large scale topology. However, since the accretion flow is fully ionized 
and freely falling, magnetic flux freezing implies that the MF-strength will 
strongly increase inwards and become in equipartition with the potential 
energy of the flow, hence terminating accretion. This impose a constrain 
on the location of the transition radius: r tr must be smaller than 10 — 20 
last-stable radii, depending on the accretion rate. 

4. Truncated disks - advective tori solution (TDAT, Hujeirat & Camenzind 2000b). 
The model relies on the results of radiative HD calculations of two-temperature 
accretion flows onto BHs without magnetic fields. These important results of 
TDAT read: 1) the disk truncate close to the last stable orbit, and forms an 
ion-dominated highly advective torus, 2) outward-oriented motions of plasmas 
between the disk and the overlying corona are formed, and which are manifested 
when MFs are included[in|. In this paper, the previous model [TH] is presented 
in details, and its astrophysical capability to reproduce the spectral energy 
distribution of the jet in M87 is confirmed. 



4 Disk structure and magnetic braking 
4.1 The outer part 

In SSDs the inward-advection of angular momentum is balanced by outward viscous 
transport. This is equivalent to equalize the terms TO and T3 in Eq. 1, while all 
other terms are neglected. 




Tl T2 T3 T4 

_^ &rB^ { ^J^r , 'l d r l 8Q | ' 1 8 ^ 8Q 
1 8i 9 88 r 2 3r ^ tur 8r cos 8 88 ° S ^ tur 88 ' 

where £ = r 2 sm8fl, r] tur = pz/ tur is the turbulent viscosity coeffecient jSH], V is the 
velocity field and B is the magnetic field (see Eq. 8). TO denotes angular momentum 
transport via advection, Tl and T2 are for magnetic extraction, and T3 and T4 
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are for viscous (micro- or macroscopic hydrodynamical turbulence) re-distribution of 
angular momentum. 

The key question here is how the constellation of these terms would look like in the 
vicinity of the black hole. 

In the outermost regions of a disk, we anticipate magnetic fields to be below equipar- 
tition. Let r tr be a transition radius (see Fig. 2), such that for r > r tr , we have the 
usual SSDs , where the ratio of magnetic to gas pressure is (3 -C 1. Here the Balbus- 
Hawley instability [2] operates on the dynamical time scale: it rapidly amplifies the 
MF and forces (3 to approach, but remains below, unity On the other hand, 

the rotational energy in SSDs exceeds the thermal energy by at least one order of 
magnitude. Therefore, the generated toroidal magnetic energy via shear can easily 
exceed thermal equipartition. Consequently, unless the imposed boundary conditions 
limit the amplification of the TMF, there is no reason to expect (3 to remain strictly 
below unity. 

4.2 The inner part 

Whether the disk is surrounding compact object or a YSO, the innermost part of the 
disk is likely to be strongly magnetized. The question which then arise is: how does 
accretion onto BHs proceed under (3 > 1 condition? 

In general, accretion via strongly magnetized and turbulence-free disks proceeds if 
the MFs are able to extract angular momentum from the rotating disk-plasma. The 
efficiency of this extraction depends, among others, on the MF-topology. For example, 
large scale or dipolar MF-topologies are considered to be appropriate for extracting 
rotational energy from the disk or even from the hole itself and power jets. In the 
present study, MFs are assumed to be of large scale topology. This is reasonable, as 
the BH-instability in combination with the Parker instability may establish a dynamo 
cycle in which the inward-advected and initially weak MFs are amplified and then 
reconnect to build up the desired MF-configuration[43j. A straightforwards conclusion 
that can be drawn is that when the matter at the outer region has completed one 
revolution around the central BH, MFs should have reached equipartition almost 
everywhere in the disk. This would eventually suppress self-generated turbulence and 
terminating angular momentum transport through turbulence-friction. 
In the region r < r tr , global conservation of the poloidal magnetic flux implies that 
Bp ~ r~ 2 . Therefore, as the innermost part of the disk contracts, there is a critical 
radius r top , below which gravitational-equipartition is maintained (i.e., Lorentz forces 
are of the same order as the gravitational forces), eventually terminating accretion. 
This configuration is unlikely, as this would imply that most of accretion inflow would 
be re-directed into outflow, with an unacceptably large mass load. On the other 
hand, if r top is sufficiently close to the last stable orbit, a change of the MF-topology 
associated with a significant loss of magnetic flux through the horizon would be 
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possible. This is a plausible event, as MFs close to the horizon are unlikely to adopt 
other than a monopole like-configuration (see Fig. 2/b). Under such circumstances, 
the generation of will be considerably reduced and the accretion of matter within 
r top proceeds unhindered. 

Around the radius r top angular momentum transport is mediated through magnetic 
braking. When combining Eq.l with the equation corresponding to the time-evolution 
of the TMF (see Eq. 8), and assuming the flow to be locally incompressible, we obtain 
a magnetic torsional wave equation which has the approximate form: 

= VlAB T . (2) 

A denotes the two-dimensional Poisson operator in spherical geometry, and Va (= 
B p /y/p) is the Alfven speed. The action of these waves is to magnetic-brakes the 
innermost part of the disk through transporting angular momentum from the disk 
to higher latitudes. Note that B p determines uniquely the speed of propagation, 
hence the efficiency-dependence of angular momentum transport on the B p — topology. 
These torsional Alfven waves (-TAWs) propagate in the vertical direction on the time 
scale: 

H 

TTAW ~ TJ- ~ Tdyn- 

This should be similar to the removal time scale of angular momentum from the disk 4 : 

r rem ~ pV T H/B P B T ~ r 3 / 2 . (3) 

Thus, T rem increases non-linearly with radius, attaining a minimum value at the in- 
ner boundary. This implies, for example, that the rate at which angular momentum 
is removed at the last stable orbit (= -Rlso) is one thousand times faster than at 
r = 100-Rlso To avoid the innermost region of the disk from running out of angular 
momentum, we require the disk to be dynamically stable. Equivalently, the rate at 
which angular momentum is removed at any radius must be equal to the rate at which 
it is advected inwards from outer layers, i.e., r a d v = T iem . This implies that the radial 
velocity — U in the disk is of the order of the Alfven speed, which is also of the order or 
even larger than the sound speed Vs, hence the terminology advection-dominated disk. 

The question which arises here is whether the disk, under these circumstances, main- 
tains its geometrical thin structure. 

The answer is hidden in the vertical momentum and energy equations. The former 

dV T , 1 dP V* VxVxB, 

— + V- W = — — + tan9-^ + L-. 4 

dt pr 89 r 4vrp w v ' 
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T rcm is obtained from the terms TO and T2 in Eq. 1 
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Taking into account that vertical transport of angular momentum is maintained 
through magnetic braking nearly without particle-advection (i.e., we neglect the sec- 
ond term in the LHS of the equation), and noting that (dP/d6)/rp is the only term 
that opposes the vertical contraction of the disk, we obtain: 



H V s 
r ~ V v 



(5) 



This implies that, as in SSDs, the thickness of the disk depends strongly on the 
temperature of the plasma, and hence on the associated heating mechanisms. On 
the other hand, since (3 > 1 in the innermost part of the disk, the dominant source 
of heating in SSDs( i.e., turbulence dissipation), is no longer efficient. In this case, 
heating is mainly due to adiabatic compression and to other non-local sources, e.g., 
radiative reflection, Comptonization and conduction of heat flux from the surrounding 
media. However, these mechanisms are unlikely to change the disk configuration 
significantly, so that the disk continues to maintain its thin geometrical structure. 
Specifically, for r < r tr , we have dH^/dr > 0. 



4.3 Formation of the ion- dominated torus 

Whether the disk model is an SSD, ADAF or ADIOS, most of the physical variable 
(e.g., density temperature, radial and angular velocities) increase with decreasing 
radius. On the other hand, the disk is said to be truncated if within a certain radius, 
say the radius of the boundary layer tbl, part of the main physical variable start 
to decrease inwards. This occurs, for example, if the central object rotates sub- 
Keplerian 5 . The Kepler-rotating and inflowing material from the disk must release 
a significant fraction of its angular momentum to be able to cross the event horizon 
of a Schwarzschild BH. This might be achieved if the MFs in the BL are sufficiently 
strong to extract angular momentum from the disk on the dynamical time scale. 
However, this requires the density to decrease considerably inwards, so that the speed 
of propagation of the TAWs largely exceeds that in the outer disk. 
We note that, unlike SSD and ADAF-models that predict a density-decrease with 
radius (i.e, p ~ r - 15 / 8 anc [ p ^ r -3/2^ reS p ec tively) , hydro- and MHD-calculations 
show, indeed, an inwards decrease of the density in the BL [IS] ,112]; [30J. 
In SSDs the pressure gradient VP gas is neglected, whereas in ADAF the pressure 
adopts the profile p ~~ r -5/ 2 . Wri en comparing the gravitational to the advective 
time scale at the inner boundary, it is easy to find that the rate of advection increases 
inwards, forcing thereby the density to be minimum at the inner boundary. As a 
consequence, the in-flowing material at r BL experiences a pressure-induced accelera- 
tion. Taking into account that the Coulomb coupling between the ions and electrons 

5 This applies for Kerr black holes with a < 1. 
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is A oc p 2 , thermal decoupling may occur, and formation of an ion-dominated torus 
is an inevitable phase in the evolution of accretion flows onto BHs, at least for sub- 
Eddington accretion. 

How does heat conduction affect the structure of the ion-torus? 
Generally, the temperature of both the ions and electrons increase inwards. The effect 
of heat conduction is to transport heat from hot into cold plasma, i.e, from inside-to- 
outside. In the absence of magnetic field and assuming the mean turbulent motion to 
be proportional to the sound speed, the innermost part of the disk will be evaporated 
by heat conduction of the ions, and therefore the torus starts to expand outwards 
When large scale MFs are taken into account, the torus contracts rather than 
expands. The reasons are: 1) heat conduction operates parallel to the MF-lines, and 
2) strong MFs act to diminish the generation and dissipation of turbulence. 
Unlike stars that heat up the surrounding media, BHs are heat-absorber (the sur- 
rounding corona is not heated from below). Therefore, the ion-dominated torus con- 
tracts and may survive in the vicinity of the last stable orbit only, where MFs are 
pre-dominantly of large scale topology. 

In addition, the location of r tr depends strongly on the accretion rate Ai. A large M. 
enhances the Coulomb coupling, reduces the effect of heat conduction and thereby 
slowing the propagation speed of the TAWs (which reduces the efficiency of the mag- 
netic braking). Thus, the volume of the ion torus shrinks and the total power injected 
into the jet-plasma decreases accordingly. 



5 The governing equations 

When modelling flow-configurations in the vicinity of a BH, it is necessary to use 
spherical geometry. Since the dynamical time scale close to the last stable orbit is 
extremely short, thermal decoupling between the electrons and the protons is possible. 
In this case, the set of equations to be solved consists of the continuity equation: 

^ + V-(pV) = 0, (6) 
the three momentum equation for the material flux: 

9M , „ /w _ T -*v , rarav.cent.rad , VxBxB 



dt 



+ V • (M (g> Vi) = -VP g + f9rav,cent,rad + V xBxB + 



where p, V — (V r ,Vg,V v ), M and P are the density, velocity vector, material flux 
M = pV, and gas pressure P(= 1Z gas p{Ti/ p x + T e /p e )). The subscripts "i" and "e" 
correspond to ions and electrons. p,\ = 1.23 and p e = 1.14. fgrav,cent,rad j g ^ Q 
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Radius (In 3xR 5CM units) 

Figure 3: A partial display of the distribution of 220 and 100 finite volume spherical 
cells in the radial and horizontal directions that covers the domain of calculations. 
Highly refined mesh distribution is adopted in the transition layer, where most of the 
physical variable undergo a strong spatial variation and where jet-launching occurs. 
The aspect ratio adopted here, i.e., the ratio of the finest volume cell to the largest 
one, is ~ 10 4 . 

force vector which includes the gravitational, centrifugal and radiative forces. Quasi- 
Newtonian gravity is used to describe the gravity of the central BH [32] • Q V1SC denotes 
the collection of second order diffusive operators. 

The evolution of the MFs is followed by solving the induction equation which reads: 

d t B = Vx{VxB-r] mag VxB). (8) 

Here, B = (B r , Bg, B v ) = (B p ,Bj>). To quench the amplification of the TMF in the 
TL, we adopt the magnetic turbulent diffusivity: ?7 ma g = aVAT-^d , where Vat = 
B'y I ^fp. In solving this equation for B , the solenoidal condition V ■ B = must be 
satisfied everywhere and for all times [T7] . 

The dynamical time scale in the innermost region of the disk is sufficiently short, 
that thermal decoupling between the electrons and ions is inevitable [SSj, [2BJ- To 
take this possibility into account, both equations describing the thermal evolution 
of the ions and electrons should be solved. In this formulation, the bulk of heat 
generated by turbulent dissipation goes into heating mainly the ions in virtue of 
their large mass compared to the electrons. The dominant cooling mechanisms for 
the ions are conduction and two-body interaction with the electrons. The electrons, 
on the other hand, are subject to various cooling processes, e.g., Bremsstrahlung, 
Comptonization, Synchrotron and conduction. Surprisingly, ion-conduction appears 
to play an important role in the innermost region of the disk, where the ions form an 
ion-dominated torus and evaporate the inner part of the disk 16 . 
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Figure 4: The horizontal distribution of the normalized angular velocity f2, the 
poloidal magnetic field components (B1,B2) = (B T ,Bg) = Bp and the TMF (Pt) 
at r = 2.5. Note the super-Keplerian rotation and the strongly enhanced strength of 
the MF-components in the TL. 



Under these conditions, the respective internal energy equations of the ions and elec- 
trons read: 

dE i 

— + V • (PV) = P*(V • V) + V • K^VTd + V - A 4 _ e , (9) 

— + V ■ (S e V) = P e (V ■ V) + V[-< ond VT c ] + A_ e _ A fl - A c - A syn (10) 

where P e ' 1 = TZ^T,/ fi u T e /fi e ), = P e >7( 7 - 1) and 7 = 5/3. V, A B , A_ c , A c , 
A syn are the turbulent dissipation rate, Bremsstrahlung cooling, Coulomb coupling be- 
tween the ions and electrons, Compton and synchrotron coolings, respectively. Here, 
we use the following express: 

3 _ _ t cn-r e ) 



Ai_ e = 5.94 x 10 niUeck 



rp3/2 



A B = 4ac/« abs p(T - E), 
k 

A c = Aan e c( -)(T e - T rad )E, 

m e c z 

where K a b s and a are the absorption and scattering coefficients. n e , rii the electron- 
and ion- number densities. The radiative temperature is defined as T rad = P 1 / 4 , 
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Radius (rn 3xR 9CH units) 

Figure 5: The distribution of the velocity field in the innermost part of the disk, 
superposed on the logarithmic-scaled ratio of the ion- to electron-temperatures (red 
color corresponds to high ratios and blue to low-ratios). 

and E is the density of the radiative energy in the zero-moment approximation of 
the radiative field. For the synchrotron cooling function, the following relation is 
adopted [34j: 

A syn = ^pAc, (11) 

where E mag , E are the magnetic and radiative energies, respectively. 

The second order operators on the RHS of Equations Q and (jTUj) are conductive 

operators and the conduction coefficients in cgs units read (Sandbaek& Leer 1994): 

n e = 7.8 x l(T 7 T e 5/2 , re, = 3.2 x l(T 8 Tf /2 

To follow the evolution of the radiative field, several assumptions have been made. 
Mainly, it is assumed that the field is isotropic and that the Flux-limited-Diffusion 
(FLD) approximation can be used to close the set of radiative moments. Thus, only 
the 0-moment equation of the radiation field is solved. Since the optical depth in the 
disk can be large in certain regions, and small in others, FLD approximation is used 
to model the radiative flux appropriately in these different regions, and to assure the 
monotonicity of the flux in the transition regions in-between. The 0-moment of the 
radiation field then reads: 

dE 

— + V • (EV) = V ■ [A FLD VE] — A B + A c + A syn , (12) 



15 



where Afld is flux limited diffusion coefficient which forces the radiative flux to adopt 
the correct form in optically thin and thick regions, i.e., 



X { V-nE if r<l, 

and provides a smooth matching in the transition regions. Here \ = p(/? a bs + cf) and 
n = VE/|VE|. 



5.1 Method of solution, initial and boundary conditions 

The solution procedure is based on using the implicit finite volume solver IRMHD3 to 
search steady-state solution for the above-mentioned 3D axi-symmetric, two-temperature, 
diffusive, radiative and MHD equations in conservative form [T7] . 
We note that the strong non-linearities governing these equation may not admit steady 
state, but strongly time-dependent or, in the best cases, quasi-stationary solutions. 
This implies that the terms describing the time- variation of the variables, i.e., d/dt, 
must be retained. On the other hand, since the dynamical time scale for accretion 
flows is shortest at the inner boundary, searching for quasi-stationary for the global 
flow-configuration on these time scales could be a computationally prohibitive process. 
To overcome this problem, we use the following spatially varying time stepping (- 
SVTS) approach, in which 

5ti = a5n, (14) 

where 5ti, 5ri are the time step size and the radial mesh increment of a finite volume 

cell at a given radius T\. The coefficient a is a constant of order unity. 

This SVTS method is appropriate for searching quasi-steady configuration of accretion 

flows or eveb for the collapse of cloud cores in molecular clouds, in which the radial 

grid points are non-linearly distributed and 5r increases with radius. If 5r is constant, 

then a must depend on the radius. In this case we suggest a = ao(r/r in ), where a is 

a constant of order unity and r in is the radius at the inner boundary. 

The main disadvantageous of SVTS method is its inability to provide time scales for 

features that possess quasi-stationary behaviour, and that might be of astrophysical 

relevance. Here we suggest to use the obtained quasi-stationary solutions as initial 

configuration and re-start the calculations using a uniform time step size. 

We note that SVTS method can be applied to explicit numerical methods as well. 

Here the global Courant number should be replaced by a local one, so that each finite 

volume cell has its own time step size. 

In the present paper, the equations are solved in non-dimensional form, using the ref- 
erence scaling variables: p = 2.5 x lCT 12 gcm _3 ,T = 5 x 10 7 K, U = V$ = j1Z g T / fr, 
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(/ii = 1.23). B = Vs^/Anp. The location of the transition layer (-TL), where the ion- 
dominated plasma is expected to rotate super-Keplerian and centrifugal-accelerated 
into jets, is shown in Fig. 2 for clarity. 

The central object is taken to be a 1O 8 M Schwarzschild BH. The domain of cal- 
culations is the quadrant: D = [0 < 9 < 7r/2] X [1 < r < 20]. Radii are given in 
units of the inner radius R in = 2.8 Rs, where Rs denotes the Schwarzschild radius 
Rs = 2GM/c 2 . Through the outer boundary and within the thickness = O.lr, 
the constant accretion rate M. = 0.01 x A^Edd is assumed. The inflow-temperature 
is taken to be T = 10~ 3 T viria i. The ion-temperature Ti is set to be equal to the 
electron temperature T e initially. A hot and tenuous corona (T = T virial , and density 
pit = 0, r, 0) = 10" 4 p(t = 0, r, 6 = 0) is set to sandwich the disk. Initially, the MF is 
set to be of large scale topology and in thermal equipartition with the disk-plasma at 
the outer boundary. Across the inner boundary, free-fall for the radial velocity and 
stress-free conditions for the angular velocity are imposed. Normal symmetry and 
anti-symmetry conditions are assumed along the equator and along the polar axis. 
The domain of integration is divided into 220 x 80 strongly stretched finite volume 
cells in the radial and vertical directions, respectively (Fig. 3). 

Special attention was given to the physical consistency of the imposed conditions that 
MFs must fulfil at the inner and outer boundaries. This is especially important in 
magnetic-diffusive plasma, as second order magnetic diffusive operators may allow 
informations to be communicated between the inner and outer boundaries on the 
dynamical time scale. Therefore, magnetic stress-free conditions have been imposed 
at both boundaries. This is equivalent to require V x B = 0. Thus, the magnetic 
diffusive operators vanish at the inner and outer boundaries, which is a reasonable 
assumption for accretion-disk inflows around BHs, in which the magnetic Reynolds 
numbers can be extremely large in certain regions and small in others. 

6 Results: formation of the super-Keplerian layer 

When an isolated Keplerian-rotating particle in the equatorial near-plane is shifted to 
higher latitudes while conserving its angular momentum, its rotation becomes super- 
Keplerian locally, and therefore it starts to move outwards. Moreover, noting that 
ttaw decreases strongly inwards, a transition layer between the disk and the hot 
tenuous corona starts to form from inside-to-outside, where the matter rotates super- 
Keplerian. Based on a previous calculations [TP] , it was found that: 1) the angular 
velocity in the TL adopts the approximate power law profile f2 ~ r - 5 / 4 ; 2) ions cool 
predominantly through fast outflows, and 3) the amplification of the generated TMF 
is equilibrated through outward advection and magnetic-quenching (i.e., we artificially 
increased the coefficient of the magnetic diffusion in the TL to allow a flux-loss of 
the TMF). In the TL, one possible solution for the radiative MHD equations is a 
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Figure 6: The distribution of the velocity field in the transition region superposed 
on equally-spaced isolines of the angular velocity Q. The strong-decrease of Q with 
radius in the equatorial region relative to its slow-decrease in the TL is obvious. 

self-similar solution in which the variables obey the following power-laws: 

fi ~ r~ 5/4 , p ~ r~ 3/2 , Ti ~ r- 1/2 , U t ~ r~ 1/4 , B P /B T ~ const. = e, (15) 

where e = H-w/r, and if\y is the width of the TL. The most prominent properties of 
this solution are: 

1. The ratio of heating to advection time scale of the ions is of order unity, i.e., 
Th/Vadv ~ 1- The ratio of heating to Coloumb-cooling time scale of the ions 
is: Th/rcoui ~ r -1 / 2 . This implies that thermal coupling between electrons and 
ions diverges with distance from their origin, and gives rise to multi-component 
jet-plasma. Furthermore, ions cool pre-dominantly through advection, i.e., adi- 
abatic expansion. 

2. The ratio of the amplification of the TMF to the flux-loss time scale is of order 
unity. The same applies to the ratio of the advection to amplification time scale. 
Thus, this implies that a significant fraction of the toroidal magnetic flux is 
advected outwards, and which thereafter acts to collimate the flow. On the other 
hand, the magnetic diffusivity in the TL scales: r/ mag ~ Hw(B p / B^)V 9 ~ er 3 / 4 . 
This 7/mag— description must break down at a certain radius r co i, where ?7 ma g 
becomes negligibly small and collimation starts to be efficient. 

3. The rate of the outward-oriented material flux in the TL reads: 

M w (r) = M w (r = r in ) + e 2 M d (r = r in )((— )^ 4 - 1), (16) 
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where .M w (r = r in ) and M.dij' = n n ) are the wind-flux rate and the disk 
accretion rate at the inner boundary. In writing Eq. 16 we have r-integrated 
the continuity equation: d Ai^/dr = Hd(pV)\g d = e 2 r ™ 3 / 4 , using H d = rsin#d, 
e = Hd/r and V = (H&/r)U. 6 d corresponds to the latitudinal angle between 
the equator and the surface of the disk is spherical geometry. 

Thus, the wind flux increases with radius, although it is a rather weak depen- 
dence. 

4. the associated angular momentum flux with the wind increases linearly with 
distance from the central BH: 

j = const, r, (17) 
which applies for r < r tr ~ 10 — 20 -Rlso 

The global energy distribution in the disk and jet region can be well represented 
by the Bernoulli number (Be). Flows with total positive energy are gravitationally 
unbound, and potentially they may propagate to infinity, whereas flows with negative 
total energy are gravitationally bound and they end their motion inside the BH. If the 
flow, however, is dissipative, energy exchange between different parts of the flow is 
possible,which gives rise to inflow-outflow configurations. Fig. 9, for example shows 
that Be is everywhere negative safe the TL, where it attains large positive values. 
We note that the outflow is sufficiently strong to shift the poloidal magnetic field 
(-PMF) lines outwards, whereas the large magnetic diffusivity prevents the formation 
of large electric currents along the equator. In the case of weak MFs (f3 < 0.1), our 
calculations reveal rather weak outflows, which is reasonable, since a low (3 yields a 
longer ttaw> which may become longer than the dynamical time scale. Moreover, low 
f3s inflows force MFs to establish a monopole like-configuration (i.e, a one-dimensional 
MF-topology which is not appropriate for producing strong TMFs). In this case, 
the rate of TMF-generation is considerably reduced. TAWs become inefficient at 
magnetic-braking the disk, reducing therefore the power required for launching jets. 
This indicates that weakly magnetized disks are in-appropriate for launching powerful 
jets, and that magnetized accretion flows are more suited for jet-production than their 
HD-counterparts. Comparing the flux of matter in the wind region to that in the 
disk, it has been found that My// Md = const. ~ e 2 . The angular momentum flux 
associated with the wind is Jyj/Jd = a(Atw/ Md) r 1 / 4 , where " db IS db constant of 
order unity. Consequently, at 500 gravitational radii almost 15% of the total accreted 
angular momentum in the disk re-appears in the wind. To clarify why the TL is 
geometrically thin, we note that the ratio: 

*™ (18) 



7"adv H A U H d 
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where r a d v is the advection time scale. The last inequality implies that angular 
momentum in the TL will be advected outwards more efficiently than being extracted 
through TAWs to higher latitudes. Furthermore, since the flow in the TL is highly 
dissipative, TMF-lines tend to close in the TL, thereby considerably reducing the 
efficiency of angular momentum transport to much higher latitudes. Also, as the flow 
in the TL rotate super-Keplerian, the normal component of the centrifugal force tend 
to shift the matter towards the equator. 

On the other hand, the corona above the TL has been found to be dynamically 
unstable. Unlike normal stars that heat up the surrounding corona from below, in 
the absence of other energy sources, black holes cannot supply the surrounding corona 
with heat, hence they start to collapse dynamically. 

To elaborate this point, let us compare the conduction time scale with the dynamical 
time scale along Bp-field at the last stable orbit of a SMBH: 

— = ^% = 4.78 x lO-Vo^fMs, 

where p±o, T^io and Aig are respectively in 10~ 10 gcm~ 3 , 10 10 K and in 10 8 Ai & units. 
This is much less than unity for most reasonable values of density and temperature 
typical for AGN-environments. In writing the above equation we have optimistically 
taken the upper limit c/y/S for the velocity, and set kq = 3.2 x 1CT 8 for the ion- 
conduction coefficient. When modifying the conduction operator to respect causality, 
we obtain r conc i/Td yn < U r /c, which is again smaller than unity. 

This agrees with our numerical calculations which rule out the possibility of outflows 
from the corona, and in particular not from the highly unstable polar region of the 
BH, where matter is neither magnetic- nor centrifugal-supported against the central 
gravity of the BH. 



7 Formation of electron-proton jets: summary 

In this paper we have presented a model for initiating jets from systems contain- 
ing BHs surrounded by accretion disks. The model is based on the following sub- 
structures: 

1. A weakly magnetized accretion disk in the outer region. Balbus-Hawley insta- 
bility amplifies the weak MFs up to thermal equipartition which, in combination 
with the Parker instability and inward-advection of matter, re-generate a PMF 
of large scale topology. This large scale PMF turns the disk at r tr into advection- 
dominated. However, since this occurs on the dynamical time scale which is 
radius-dependent, r tr is basically the outer boundary, or alternatively is close to 
the last stable orbit. The former possibility is observationally inconsistent, as 
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this would force the disk to truncate at large radii and terminate accretion. The 
later possibility is plausible: assume that the flow at r tr is freely falling, and has 
the vertical width of H d /r = 1/30. Using p(r,9 = 0) = p(r ,9 = 0)(r B L/V) 3/2 
and B p (r, 8 = 0) = B p (ro, 9 = 0)(rBL/^) 2 , one can easily verify that MFs would 
terminate accretion if r tr > 16rLso, where tlso is the last stable radius. Here 
we have assumed that whatever is the strength of the MF, the magnetized inflow 
interior to r L so is gravitationally bound and will be accreted, together with the 
MF, into the hole. In this region, we anticipate large scale MFs to change into 
a monopole like-configuration. Thus, an advective disk threaded by large scale 
MFs must be located inside 10 — 20rLso, depending on the accretion rate. 

2. In the region r top < r < r tr , TAWs are the dominant angular momentum 
carrier. The transport proceeds vertically and on the dynamical time scale, 
which requires the disk, again, to be advection-dominated for stably supplying 
the jet with angular momentum. The angular momentum associated with TAWs 
is deposited in the turbulent-diffusive TL between the disk and the corona, 
where vast of the shear-generated TMFs reconnect and heat the ions up to the 
virial temperature. Besides the fact that the large magnetic-diffusivity in the 
TL lengthens the time scale of TAWs to cross the TL compared to the local 
dynamical time scale (which is necessary for maintaining a geometrical thin 
TL), reconnection of the TMF-lines terminates magnetic braking and yields a 
trapping of the angular momentum. Consequently, the virial-heated ions are 
then forced to rotate super-Keplerian and centrifugal-accelerate outwards. 
Another argument in favor of the thin geometrical structure of the TL is that 
the corona is thermally unstable. The reasons are: 1) absence of heating from 
below, 2) centrifugal forces are weak, and 3) Lorentz forces exerted on the 
plasma in the polar region of the BH are negligibly small. 

3. The plasma in the TL is tenuous, highly advective, two-temperature and ion- 
dominated, which gives rise to the formation of electron-proton jets. The shear- 
generated TMFs in the TL are in thermal equipartition with the ions and in 
super-equipartition with the thermal energy of the electrons. Such a strong 
TMF in the TL is necessary for explaining the observed excess of radio lumi- 
nosities such as in M87 and GRO 1655-40. 

4. Since the pressure-gradient vanishes at the horizon, a rarifaction wave starts 
to move from inside-to-outside, leaving the matter behind its front free-falling. 
The resulting density profile has a global maximum at r B L(<< r tr ), and a min- 
imum at the inner boundary r in , where the dynamical time scale is shortest, 
and where the strength of Coulomb coupling between the ions and electrons 
is weakest. Consequently, the disk truncates at r B L and an highly advective 
ion-dominated torus emerges in the BL. The MF in the ion-torus change topo- 
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logically from large scale into monopole. Such a change is associated with a 
significant loss of magnetic flux in the hole, which is necessary for preventing 
magnetic termination of accretion. 

5. When does the outflow start to collimate into jet? 

A significant part of the shear-generated TMFs reconnect in the TL and heat 
the ion preferentially. The other part is advected outwards with the plasma on 
the dynamical time scale. Inspection of Equation 1 shows that the ratio of the 
viscous to the advective time scale of angular momentum in the TL reads: 

^ « (19) 

7"adv ^tur 

In writing this approximation we have inserted V^" = B^j ^fp = V v , which is 
applicable for the plasma in the TL. Therefore, given the critical radius r co i(= 
ftur/VA) at which both viscous and advective time scale are equal, collimation 
will start for r > r co \. This is easy to fulfil, as Va increases with distance from 
the BH, or remain constant, whereas z/ tur must decrease. This gives rise to 
extremely large ratio of r visc /r adv . Therefore, beyond r co \, particle-trajectories 
of motion are dictated by MFs. 

Finally, the model is similar to the Blandford & Znajeck (1977) scenario, as vertical 
transport of angular momentum by TAWs is mediated almost without advection of 
matter. This gives rise to high collimated and extremely low mass-loaded jets. The 
global configuration of the jet-disk connection presented in this paper enhances the 
role of the central engine in powering the jet [21]. A spinning BH will likely feed 
the jet with extra rotational energy through the frame dragging effect, whereas a 
non-rotating hole may have the opposite effect: extraction of energy from the plasma 
in the TL, hence a less powerful jet. 

Additionally, the model may provide explanations for several other important issues 
related to the formation of jets and disks: 1) the Ka emission lines observed in Cyg 
X-l which suggest that the disk-plasma in the vicinity of the last stable orbit is cold 
and rotates Keplerian, 2) The origin of the super-Keplerian motions observed in the 
Galactic black holes XTE J1650-500 and MCG-6-30-15 which suggest that the 
central objects are spinning BHs, and 3) the origin of the extremely low mass-loaded 
jets observed in several AGNs and /i-quasars. 
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Figure 7: The distribution of 30 weak poloidal magnetic field lines around a non- 
rotating BH after 10 orbital periods at the inner boundary (black lines). The TMF- 
distribntion is shown as well (bine color corresponds to lar^e nositive values, and 
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Figure 8: 30 equally-spaced iso lines of the poloidal-component Bp ( black lines) and 
the density-distribution (yellow color corresponds to very high density-values, blue to 
middle and violet to low values). 
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Figure 9: A snap-shot of the distribution of the Bernoulli number in two-dimensional 
quasi- stationary calculations. The decrease from large to low positive values is rep- 
resented via yellow, green and red colors. The blue color corresponds to negative 
values. The figure shows two ejected gravitationally unbound blobs of large positive 
energies in the TL. 
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Figure 10: The SED of a disk-jet model of the elliptical galaxy M87. Similar to 
Biretta 2002 and Di Matteo 2003, an accretion rate of 5 x 10~ 4 Aldington and T = 
5 x 10 6 K is set to enter the domain of calculations through the outer boundary which 
is located at 150 Schwarzschild radii from the central SMBH. The vertical scale height 
of the disk at the outer boundary is taken to be H rs 0.1i? out . Additionally, a hot 
tenuous corona is set to sandwich the optically thin disk. The calculated profiles (solid 
lines) are superposed on the observational data (asterisks). The line 07 corresponds 
to a model in which the PMF is set to be in equipartition with the thermal energy, 
whereas TMF=0. The line 08 is similar to 07, but the TMF is allowed to develop and 
reach values beyond equipartition with respect to the thermal energy of the electron 
in the TL. The above SED has been obtained by solving the radiative transfer in 4- 
dimensions, taking into account the Kompaneets equation for consistently modelling 
Comptonization (see [H^EO] f° r further details). 
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